Evidence for deconfined quantum criticality 
in a two-dimensional Heisenberg model with four-spin interactions 
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Using ground-state projector quantum Monte Carlo simulations in the valence bond basis, it is 
demonstrated that non-frustrating four-spin interactions can destroy the Neel order of the two- 
dimensional S = 1/2 Heisenberg antiferromagnet and drive it into a valence-bond solid (VBS) 
phase. Results for spin and dimer correlations are consistent with a single continuous transition, 
and all data exhibit finite-size scaling with a single set of exponents; z = l,v = 0.78 ± 0.03, and 
77 = 0.26 ± 0.03. The unusually large 77 and an emergent (7(1) symmetry, detected using VBS order 
parameter histograms, provide strong evidence for a deconfined quantum critical point. 

PACS numbers: 75.10.-b, 75.10.Jm, 75.40.Mg, 75.40.Cx 



Since the discovery in 1986 of high-Tc superconductiv- 
ity in layered cuprates, quantum phase transitions in two- 
dimensional (2D) antiferromagnets have formed a central 
topic in condensed matter physics [J, Q • While supercon- 
ductivity is induced in the Cu02 planes of the cuprates 
by doping with charge carriers, other mechanisms for de- 
stroying the Neel order and stabilizing different ground 
states have also been intensely investigated theoretically. 
Considerable efforts have been devoted to possible spin 
liquid ("RVB" and valence-bond solid (VBS) states 
driven by magnetic frustration 0, H, Q • This work has 
been partially motivated by the hope that an understand- 
ing of generic features of quantum phase transitions in 
2D antiferromagnets could shed light also on the mecha- 
nisms at work in the cuprates Quantum fluctuation 
driven phase transitions are also of broader relevance in 
the context of strongly correlated systems Q. 

A quantum phase transition occurs as a function of 
some parameter at temperature T = and corresponds 
to a T > transition in an effective classical system with 
an imaginary-time dimension — the path integral Q . The 
standard Landau-Ginzburg- Wilson framework for critical 
phenomena should thus be applicable, with the dimen- 
sionality rf ^ d + z , where the dynamic exponent z de- 
pends on the way space and time correlations are related. 
In the paradigm prevailing until recently, the "Landau 
rules" for the nature of the transition — continuous or 
first-order — were also assumed to remain valid for quan- 
tum phase transitions. A direct transition between two 
ordered phases should thus be generically first-order if 
two different symmetries are broken. This notion has 
recently been challenged by Senthil et ai, who argued 
that quantum phase transitions separating two ordered 
phases can be generically continuous, even when differ- 
ent symmetries are broken jlC|. This theory of "decon- 
fined" quantum critical points was first developed for 
the transition between an antiferromagnetic (AF) and a 
valence-bond-solid (VBS) state. Both these states have 
confined S = 1 excitations — gapless magnons and gapped 
"triplons", respectively. The critical point is character- 



ized by deconfined S" = 1/2 spinous coupled to an emer- 
gent U{1) gauge field [l3|- In 2D the deconfined state is 
unstable and exists only at a point separating the two or- 
dered phases. The AF and VBS order parameters arise 
as a consequence of spinon confinement. In this Let- 
ter, quantum Monte Carlo (QMC) results are presented 
which support this theory. 

Preceding the theory of deconfined quantum critical 
points, continuous transitions between two ordered quan- 
tum states had been suggested based on numerical sim- 



ulations [ill, II4I. However, in more detailed studies fol- 



lowing the theoretical developments it has proved difh- 
cult to confirm their existence. Instead, many studies 
have pointed to weakly first-order AF-VBS transitions 
H M, E H m or other scenarios inconsistent with 
deconfined quantum criticality [3]. To date, large-scale 
QMC studies of potential deconfined quantum critical 
points have focused on spin (or hard-core bosoriic) models 
with spin-anisotropic interactions 1^, 14, 1^ lg|. Frus- 
trated SU (2) (Heisenberg) symmetric interactions, which 
cannot be studied using QMC simulations due to the in- 
famous " sign problem" , have been considered in exact 
diagonalization studies [l9|. Because of the limitations 
to very small lattices, it has not been possible to study 
phase transitions in detail, however. In fact, not even 
the nature of the VBS state has been completely settled 
in basic models such as the J1-J2 frustrated Heisenberg 
model [20|]. 

Here it will be shown that the AF order of the square- 
lattice Heisenberg model can be destroyed also by non- 
frustrated isotropic interactions accessible to QMC sim- 
ulations. The following Hamiltonian will be discussed: 



i7 = J^S,.S, -Q J2 (S,-S,-i)(Sfe-S, 

{ij) (ijkl) 



(1) 



where (ij) denotes nearest-neighbor sites and (jjfc?) refers 
to the corners of a plaquette, such that ij and kl form two 
parallel adjacent horizontal or vertical links. This inter- 
action contains a subset of the four-site ring-exchange, 
and with Q > there is no QMC sign problem. Note 
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FIG. 1: (Color online) Finite size scaling of the squared spin 
(M) and dimer (_D) order parameters at J/Q = and 0.1. 
The curves are cubic fits. Statistical errors are much smaller 
than the symbols. 

that the purpose here is not to model any specific mate- 
rial, but simply to construct a model system in which an 
AF-VBS transition can be investigated. It will be shown 
below that the ground state of the J-Q model has AF 
order for J/Q > 0.04 and VBS order for J/Q < 0.04. 

To study the ground state of the hamiltonian ([T]), an 
approximation-free projector technique in the valence 
bond basis 2l| is employed which is ideally suited for 
multi-spin interactions formed out of singlet projection 
operators (Sj • Sj — j). Here L x L lattices with L up 
to 32 are considered. Larger systems may be reachable 
using loop algorithms in the standard basis, which 
have been used for U{1) models with four-site interac- 
tions [23, [i^- The valence bond basis has its advan- 
tages, however, including an improved estimator for the 
singlet-triplet gap. 

Results will be presented for spin-spin (s) and dimer- 
dimer (d) correlation functions, 



Cs{r) - (S(O).S(r)), 

C4r) = ([S(0).S(i)][S(r).S(r- 



x)]). 



(2) 
(3) 



where x denotes a lattice unit vector in the x direction. 
The AF order parameter is the staggered magnetization, 
the square of which is calculated; 



(4) 



The VBS state can have either columnar or plaquette 
order, both of which break Z4 symmetry. An impor- 
tant aspect of the theory is that these order parameters 
should both exhibit divergent fluctuations at the decon- 
fined critical point. Only at some length-scale diverging 
as a power of the correlation length should one of them 
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FIG. 2: (Color online) Finite size scaling of the singlet-triplet 
excitation gap multiplied by the system size L. The behavior 
at J/Q = 0.04 corresponds to z = 1. 

be singled out • This is analogous to the irrelevance of 
Z4 anisotropy in the 3D XY model [2^ and corresponds 
directly to the predicted emergent U{1) symmetry. The 
q = (tt, 0) dimer order parameter, 



(5) 



is divergent for both columnar and plaquette VBS order 
and will be studied here. 

Extrapolations of the AF and VBS order parameters, 
shown in Fig. [l] demonstrate that there is long-range 
VBS order but no AF order at maximal four-spin inter- 
action; J/Q = (note that there are still two-site inter- 
actions present when J — 0; simulations for J < are 
sign problematic). Also shown are results at J/Q — 0.1, 
where the situation is the reverse; there is AF order 
but the VBS order vanishes. Thus there is an AF-VBS 
transition somewhere in the range < J/Q < 0.1, or 
there could be a region of AF/VBS coexistence (which 
would be analogous to a supersohd state). The nature 
of the VBS order — columnar or plaquette — is not clear 
from these calculations. However, simulations of open- 
boundary rectangular lattices, in which a unique colum- 
nar or plaquette pattern can be stabilized [12i] . indicate 
that columnar order is preferred. The extrapolated VBS 
correlation at J/Q = is sa 0.0024. 

The deconfined theory has dynamic exponent z — 1 
|l3] ■ This exponent can be directly accessed through the 
finite-size scaling of the singlet-triplet gap; A ^ L'". 
To demonstrate consistency with z = 1, the scaling of 
LA is shown in Fig. [5] for J/Q — and 0.1, as well as 
for J/Q = 0.04 which will be shown below to be close 
to criticality. Here LA extrapolates to a non-zero value, 
supporting z = 1, and at J/Q = and 0.1 the behaviors 
are what would be expected off criticality. The inset of 
Fig. [2] shows an infinite-size extrapolation of the gap at 
J/Q = 0, giving A/Q « 0.07. 

Correlation lengths and for spins and dimers are 
defined in the standard way as the square-roots of the 
second moments of the correlation functions ^ and 
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FIG. 3: (Color online) Scaling with u = 0.78 and = 0.04 of 
the correlation lengths (a) and the spin Binder ratio (b). 

Also useful is the Binder cumulant, defined for the spin as 
Qs = (M^)/(A/^)^. Finite-size scaling of these quantities 
is used to extract the critical coupling and the correla- 
tion length exponent ly. To achieve good data collapse, 
a subleading correction is also included. With g = J/Q, 
the scaling ansatz is, 



Aig,L) ^ L^il + aL--)f[ig - g,)L'/'^] 



(6) 



where A = ^s^S,d, or g^, and k = 1 for ^siCd ^-nd for 
Qg. As seen in Fig. [31 these quantities can be scaled using 
gc = 0.040 ± 0.003 and a common i/ = 0.78 ± 0.03. In all 
cases, the subleading exponent a; « 2, and the scaling is 
nearly as good if w = 2 is fixed throughout. Interestingly, 
the best prefactor a is then almost equal for and ^d, 
a « 8, but this may be coincidental. 

Next, the correlation functions Cs,d{^) at the longest 
lattice distance, r = (L/2,L/2), are analyzed to extract 
the correlation function exponent 77. The expected scal- 
ing is as in Eq. ^ with k = — (1 + r/). Now g^ and 1/ 
are kept fixed at the values determined above. As shown 
in Fig. [31 a single exponent describes both the spin and 
dimer data, and in this case a subleading correction is 
not needed (a = 0). The exponent, rj = 0.26 ± 0.03, 
is unusually large. In the 3D 0(3) universafity class, 
describing transitions between the AF and a featureless 
gapped state 0, [11], rj ~ 0.04. A larger rj for deconfined 
quantum criticality was argued for on physical grounds 
by Senthil et al. [10|. The universality class was argued 
to be that of the hedgehog suppressed 0(3) transition, 
for which 13 /ly = (1 + 77)72 0.80 ± 0.05 was obtained 
in simulations of a classical model in [2^. This is larger 
than P/v = 0.63±0.02 found here, but on the other hand 
smaller lattices were used in [2^ and there may also be 
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FIG. 4: Long-distance spin and dimer correlations scaled us- 
ing V = 0.78, T] = 0.26, and = 0.04. 

issues with how hedgehogs were suppressed. The direct 
study of an actual AF-VBS transition presented above 
can thus be expected to be more reliable. 

It is also interesting to study the probability distribu- 
tion P{Dx, Dy) of the dimer order parameter. In the VBS 
phase, one would expect this distribution to rcfiect the Z4 
symmetry, i.e., for a columnar VBS there should be peaks 
at Dx — 0, Dy — ±D and = ±1?, Dy ~ (whereas a 
plaquette state would give rise to peaks rotated by 45°). 
It should be noted, however, that P{Dx,Dy) is a basis 
dependent function. In the valence bond simulations (2ll | 
the order parameters used to construct P{D.j;,Dy) are 
matrix elements (with e = x,y). 



Dr^ 



^b|^E,.S(r)-S(r + e)(~l)'--|v['.) 



(7) 



where |5'a), |^b) are valence bond states generated by 
operating with a high power on trial state (stochas- 
tically sampling valence bond evolutions). Although 
P{DxT Dy) is not a physically measurable quantity, any 
symmetry detected in it should refiect an underlying sym- 
metry of the projected state. Fig. [5l shows a color-coded 
P{Dx,Dy) histogram generated at J/Q = 0. The ex- 
pected Z4 symmetry of the VBS is not seen; instead the 
histogram is ring shaped, which indicates a U (1) symmet- 
ric order parameter. Such an emergent U{1) symmetry 
is in fact predicted by the deconfined theory in the 
VBS phase below a length scale A which diverges faster 
than the VBS correlation length; A ~ with a > 1. 
Thus, inside the VBS phase, if the system length L <C A 
one should expect to find an C/(l) symmetric order pa- 
rameter, with the Z4 becoming relevant only for larger 
sizes (and then seen as four peaks emerging in the his- 
togram). Here, apparently, even at J/Q = the system 
is close enough to the critical point for the system length 
{L — 32) to be less than A and, hence, Z4 to be irrele- 
vant. Recalling that the VBS gap is small, A/Q « 0.07, 
and that A ~ A~°, this seems reasonable. On 

moving closer to the critical point, P{Dx,Dy) smoothly 
evolves into a single broad peak centered at (0,0), as 
is expected for a continuous transition. Note that the 
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FIG. 5; (Color online) Histogram of the dimer order param- 
eter for an L = 32 system at J/Q = 0. The ring shape 
demonstrates an emergent (7(1) symmetry, i.e., irrelevance of 
the anisotropy of the VBS order parameter. 

finite-size extrapolation of the order parameter in Fig. [T] 
is not sensitive to the nature of the VBS state — plaquette 
or columnar — and should give the correct magnitude of 
the order parameter even though no Z/^ features are yet 
seen in the histogram for these system sizes. 

The above analysis points consistently to a deconfined 
quantum critical point as the most likely scenario for the 
AF-VBS transition in the J-Q model. One set of ex- 
ponents describes both spin and dimer correlations, the 
value of T] is unusually large, and there is an emergent 
U(X) symmetry in the VBS order parameter. In prin- 
ciple one cannot rule out a weakly first-order transition 
on the basis of finite-size data. However, although the 
lattice sizes used in this work are not extremely large, 
it would be hard to explain why a first-order transition 
should lead to the kind of scaling observed. A narrow 
region of AF/VBS coexistence is also unlikely, as there 
would then be two transitions and there is no reason to 
expect the spin and dimer critical exponents to be the 
same (in particular, the unusually large j]). It is difficult 
to say anything more quantitative regarding a possible 
first-order transition or coexistence based on the calcu- 
lations presented here. 

An emergent C/(l) symmetry may also explain why 
it has been so difficult to determine the nature of the 
VBS state in the J1-J2 Heisenberg model [23|. Even if 
the transition would be weakly first-order in this case 
27[, an emergent UiV) symmetry could still affect 
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small lattices, thus making it difficult to distinguish be- 
tween columnar and plaquette VBS patterns. Emergent 
?7(1) symmetry may be more common than deconfined 
quantum criticality and could hence affect many models 
with VBS states. The high density of low-lying singlets 
associated with C/(l) symmetry may also affect exact di- 
agonalization studies [a] of level spectra. 
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